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ABSTRACT 


This thesis presents a propagation study and antenna design for the 
Micro-Remotely Powered Vehicle (MRPV). A propagation loss study was 
conducted to determine the attenuation of various building walls and to select 
an optimum frequency band for antenna design. An approximate ray tracing 
model for loss determination was developed and programmed in MATLAB. The 
computed losses from the model are presented for comparison with measured 
results. One possible antenna configuration for the MRPV, a circumferential 
slot, is analyzed and design parameters varied to obtain optimum antenna gain. 
A prototype of the slot antenna was developed. The antenna patterns, 


efficiency, and Voltage Standing Wave Ratio (VSWR) are presented. 
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l. INTRODUCTION 


A. BACKGROUND 

1. Beamed Microwave Power 

The concept of beamed microwave power transmission was first 
introduced in the 1960’s. Under U.S. Air Force funding, the Raytheon Company 
led the pioneering research efforts which resulted in an S-band (2.45 GHz) 
power beamed helicopter platform in 1964 [Ref. 1]. The demonstration of a 
beam riding helicopter soon followed in 1969. The beam riding helicopter 
guided itself using a microwave beam as a position reference for roll, pitch, yaw, 
and x and y translation. As a means of transferring energy, the main advantage 
of beamed microwave power transmission is in placement of the power source 
on the ground. A ground based microwave power source allows for the design 
of a very small platform with unlimited duration. Due to the large mass to power 
ratio required for batteries and internal combustion devices, it is not possible to 
design very small airborne platforms for even short ranges given current 
technology. A diagram of the principal elements for beamed power 
transmission is shown in Figure 1. 

The beamed power transmission consists of a dc to microwave power 
conversion, transmitting antenna, free space transmission with possible 
blockage, and the microwave power reception and conversion to de. The 
theoretical maximum efficiencies for each element of transmission are typically: 
dc to microwave power conversion - 90%, antenna efficiency - 97%, reception 


and conversion to dc - 92%. Therefore, a maximum overail dc to dc efficiency of 





| dc to Microwave Power 
Microwave Antenna Reception and 


Power Conversion to dc 
Conversion 





Figure 1. Beamed power transmission 


76% is possible. A certified experimental overall dc to dc efficiency of 54% has 
been achieved in the laboratory [Ref. 2]. 

The beamed power is limited in transmission by the loss resulting from 
the spreading of the electromagnetic wave front as it propagates with distance 
and the loss due to any blockage encountered. The spreading loss is a function 
of distance, and beyond that, can only be improved by increasing the directivity 
of the antenna field pattern. The blockage loss depends on the electrical 
properties and physical characteristics of the obstruction. Therefore, the 
transmitted power must be determined from the required receive power and 
losses expected during the propagation of the beamed microwave power. 

The reception and conversion from beamed microwave power to dc is 


performed by a device known as the “rectenna.” The rectenna, as its name 























suggests, consists of an antenna and a rectifier. In its simplest form, the device 
can be a half-wave dipole with impedance matching circuitry, terminating into a 
rectifying diode. 

2. Micro-Remotely Powered Vehicle (MRPV) 

The MRPV is a microwave power beamed surveillance helicopter 
currently under development with the Advanced Research Projects Agency 
(ARPA). The system requirements for the MRPV are as follows: 

1) Provide immediate, continuous, video surveillance of the inside of any 

desired building or volume to a nearby unit. 

2) Remain inexpensive and be covert. 

3) Maintain Department of Defense (DOD), Department of Justice (DOJ), 

and commercial applications. 

The MRPV is a very small (~2 cm2), quiet, rotary wing aircraft with 
unlimited mission duration using a video camera for both real-time platform 
control and telemetry. The MRPV will be able to fly out at ranges greater than 
50 m and monitor a scene with unlimited mission duration. The small size will 
also make the MRPV difficult to detect or destroy. The applications for MRPV 
will include: indoor building surveillance, monitoring of disputed borders, 
tactical military surveillance, and commercial remote monitoring. The MRPV 


geometry is shown in Figure 2. 





— Se 5 com Cs 
Main Rotor 
ne” Se \ 
—_— 


RF Antenna Area 


Camera Lens 


Volume = 6 cm? /<¢——— 2cm —_— RF Antenna Area = 3 cm@ 


Weight = 35 g Lens Diameter = 0.4 cm 





Figure 2. MRPV geometry 


B. OBJECTIVE 

The objectives of this research are: 1) to investigate the propagation 
losses of typical building walls, and 2) to design, build and test a prototype 
antenna for reception of beamed microwave power for the MRPV. 

The propagation loss study will quantify the losses of a microwave power 
beam through typical building walls, the results of which will be used to select a 
system operating frequency and transmit power level. The type and location of 
the building walls used in this study are presented with the actual loss 
measurements. The loss measurements were taken using a HP 8350B signal 
generator and a HP 9562A spectrum analyzer. There are two primary 


frequencies of interest which arise from requirements to use existing off the 











shelf hardware for the transmitter. These are 2.5 GHz (scientific and 
commercial band) and 5.5 GHz (satellite uplink frequency). The loss 
measurements were taken over a frequency band from 2 to 6 GHz. 

An approximate model is presented for the propagation of a spherical 
electromagnetic wave front through a wall. This model is based on the method 
of images and geometrical optics. The propagation of locally plane sections of 
a spherical wave front are used to derive an equivalent free space source (i.e., 
with the wall removed) to account for the phase change and loss differences 
along the spherical wave front. Because of the locally plane assumption, the 
Fresnel transmission coefficients for a plane electromagnetic wave front 
passing through a lossy dielectric can be used. The transmission through a 
wall includes both the direct and floor reflected components of the wave. An 
equivalent source is derived for the direct waves and an equivalent image 
source is derived for the reflected waves. The complex values of the direct and 
reflected fields are computed at points representing the antenna aperture. The 
resulting array of field points is then integrated to find the total field received by 
the antenna. 

The results of the MRPV antenna design are also presented. The 
primary antenna configuration of interest for a vertically polarized electric field is 
the circumferential slot. Calculations of far-field patterns and surface currents 
are simulated using the PATCH computer code which uses the method-of- 
moments (MM) solution to the electric field integral equation. Based on the 
PATCH code results, a prototype antenna was designed and built. The far-field 
patterns and efficiency were measured in the Naval Postgraduate School 


(NPS) anechoic chamber for a scaled model of the antenna. An approximate 


antenna efficiency of 67.5% was measured. Using known values for power 
conversion efficiency, a possible rectenna configuration is shown with an 
overall microwave to de power conversion efficiency of 61%. 

The remainder of this thesis consists of four chapters. Chapter Il 
discusses the first part of the propagation loss study. The measurement system 
for determining the propagation loss of various walls is presented, followed by a 
description of the walls tested and loss measurements. Chapter Ill continues 
the propagation loss study presenting the approximate ray tracing model that is 
derived for comparison with measured results. The model is based on the 
method of images and geometrical optics and approximates the propagation of 
a spherical wave front through a lossy dielectric boundary. Chapter IV presents 
the antenna design. A simulation of the prototype antenna is performed to 
determine the optimum feed point. Chapter V presents the results of 
measurements performed on the antenna. This includes the antenna 
parameters and the far-field measurements taken in the anechoic chamber. 
Chapter VI provides conclusions and any recommendations for follow on 


research in this area are also noted. 

















ll. PROPAGATION LOSS MEASUREMENTS 


A. MEASUREMENT SYSTEM 

The wall loss is estimated by measuring the power received from a 
transmitter when the wall is present and comparing it to the received power 
without the wall. The propagation loss measurements were taken using the 
system depicted in Figure 3. The system consists of transmitting and receiving 
stations each with parabolic reflector antennas. The transmit signal is provided 
by a HP 8350B signal generator operating in the frequency sweep mode. The 
signal is received using a HP 8562A spectrum analyzer utilizing the trace hold 
feature. Used in this configuration, the transmitter sweeps through a specified 
frequency range with the receiver storing a trace of the resultant signal strength. 
Baseline free space measurements are taken to serve as a reference. The free 


space reference data was measured at separation distances of 10 and 15 ft. 


Meter 
Frequency eo Plotter 


Counter 
| Signal Generator | Signal Generator 


= |Z 


Transmitting Receiving 
Station Station 





Figure 3. Propagation loss measurement system 





1. Transmitting and Receiving Stations 
Block diagrams of the transmitting and receiving stations are given in 


Figure 4. 
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Figure 4. Measurement system block diagram 


As shown, the transmitter station includes a HP 8350B signal generator, 
HP 5350B frequency counter, and a HP 436A power meter. The readings for 
frequency and power are obtained using 20 dB directional couplers as shown. 
The transmitting antenna is connected to the system using approximately 6 ft of 
50 ohm coaxial cable. The receiving station consists of an HP 8562A spectrum 
analyzer and HP 7475A pen plotter connected together on a standard HP-IB 


interface. The receiving antenna is connected to the system using 























approximately 7 ft of 50 ohm coaxial cable. The transmit and receive antennas 
are both parabolic reflectors with log periodic feeds. 

The free space reference data is stored in the HP 8562A spectrum 
analyzer for subtraction with the blockage measurement. It is stored using the 
A-B-->A math function of the spectrum analyzer. The spectrum analyzer is also 
set in a maximum hold state for the A-B-->A trace. The trace records the 
measured propagation loss as the transmitter sweeps across the set range of 
frequencies. 

2. Transmit and Receive Antennas 

Identical transmit and receive antennas were assembled using American 
Electronic Laboratories (AEL) model APN102C parabolic reflectors with AEL 
model APII-101B log periodic feeds. The antennas were designed to be 
portable using modified Tektronix K212 instrument carts. The design and set- 
up of the portable parabolic reflector antennas is shown in Figure 5. In this 
figure, the antenna dimensions and electric field orientation with respect to the 


log periodic feed are also provided. 


Tektronix K212 


Instrument Cart 





Figure 5. Parabolic reflector antenna with log periodic feed 





3. Equipment Listing 
A listing of all equipment used for both the transmitting and receiving 


stations is provided in Table 1. 






(Serial # [Description 
oo 
P 83508 0.005 - 12 GHz CW, -2 to +18 dbm 
P 53508 
PA36A | T8OSA0ZEA 
API 1078 

A|-_____|S00hm 0st 18GHz 
G 27470 
cr 
P7790 
P7790 
hs nr 
eis ane 
HP O562A 
P7475 
API 107 
RG 9670 
AG O70 
HPIB___[-____[analyzertoplotierinterface 


Table 1. Equipment listing 
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B. MEASUREMENTS 

The system described in Chapter Il, Section A, was used to measure 
propagation loss of a microwave power beam through several building wails. 
The measurements were taken using carefully chosen reference marks on the 
antennas for equipment alignment. With a standardized alignment of the 
antenna platforms and procedures, accurate comparisons can be made 


between measurements taken with and without the wall. The loss 
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measurements were collected over a frequency band of 2 to 6 GHz. The test 
measurement configuration is depicted in Figure 6. The distances d; and do 
are those from the wall to the transmitter and receiver, respectively. The wall 
has thickness L, and material constants e, (dielectric constant), Lo 
(permeability), and tan 5 (loss tangent). The distances from each antenna to the 
wall were set equal (dj=d2) yielding a total distance dy+L+do. Transmit and 
receive distances of 10 and 15 feet were used. These distances were the 


maximum possible given the confinement of building halls and rooms. 


Eo,Ho Er, Uy, tan d 


Ap IN 


o> | — 
CLIO PP ALP DALE EEE AE ED EL LLLLLAA MAA MMM aRAAAAAAAA ALA AM hh 


dy L do 





Figure 6. Test set-up 


1. Building Walls 

The walls measured in this test are located on the campus at the Naval 
Postgraduate School. A brief description of each wall is given in Table 2. The 
exact composition of the walls is unknown. The descriptions in Table 2 are 


based solely on visual inspection. 
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Building Wall Material Type 


Front Entrance 

1.75 " Double Doors Spanagel Hall, 
Room 429 
w/ Chalkboard cement Room 408 
Wall Room 421 

glass windows Roof exit doors 


Table 2. Typical building walls 











—_h 


2. Measured Data 
The measured losses for the materials listed in Table 2 are provided in 
Figure 7. The measured data for each wall was plotted using the HP 7475A 


pen plotter and the raw data is included in Appendix A. 


Measured Propagation Loss —#— A- .25" Glass 


5 —O— B- 1.75" Wood 
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10.75" Cement 
Wall 


Loss (dB) 


1.75" Metal 
Freq (GHz) Covered Doors 





Figure 7. Loss measurements of typical building walls 
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ill. PROPAGATION LOSS MODELING 


A. BACKGROUND 

An approximate model was developed for the propagation of a spherical 
electromagnetic wave through a building wall of known material properties. 
This model is based on the method of images and geometrical optics. The 
propagation of locally plane sections of a spherical wave front is used to 
determine the signal at the receive antenna that includes the phase and 
amplitude changes introduced by the wall. Based on the received signal, an 
equivalent free space source can be defined and the wall removed. A similar 
process is used to determine a free space image source located below the floor. 
Thus, these two sources radiate in free space and their superposition yields a 
solution that includes both the direct and reflected transmit waves. The complex 
values of the direct and reflected fields are computed at points representing the 
antenna aperture. The resulting array of field points is then integrated to find 
the total received field of the antenna. This approximate model is shown to be 
in good agreement with measured results. 

1. Spherical and Planar Electromagnetic Waves 

A spherical electromagnetic wave front emanates from an ideal point 
source. The ideal point source refers to an element of uniform amplitude 
current which is of infinitesimal length. The electric field vector of all realizable 
antennas contains r, Q, and 6 field components. The r component of E cannot 
be assumed negligible when operating in the near-field region of the antenna. 
In the far field the radiated electromagnetic wave is spherical. Furthermore, 


only the @ and 6 components are significant and they vary as Yr. Locally the 
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spherical wave can be considered to be a plane wave, when r is large as 


shown in Figure 8. 


Spherical wave with Plane wave 
large radius of curvature approximation 





Figure 8. Local plane wave approximation 


2. Plane Waves at Oblique Incidence on an Interface 

The transmission and reflection of a plane wave at oblique incidence to a 
boundary is shown in Figure 9. In this derivation the boundary is assumed to be 
between free space and a lossy dielectric. The medium is characterized by ¢,, 


u,, and tan 6. 


Er, Ur, tan 6 


Re 


Ot 





region 1 region 2 





Figure 9. Plane wave at oblique incidence - single interface 
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The perpendicularly polarized wave, consisting of E in the plane of 


incidence and H parallel to the boundary, is impinging on a boundary at an 


angle of incidence 6;. The boundary condition requires that [Ref. 3] 


B, sin®; = By sin®, = Bo sin®, (3.1) 


or, 8, = 9;,, which is Snell's Law. The phase constants are given by Bo =2n/Xdo 
and Bo =2n/Ap. The wavelengths in free space and the medium are denoted 


Xo and Az, respectively. Knowing that the intrinsic impedance of the dielectric 


regions are given by [Ref. 1] 





No = =n (region 1) (3.2) 
0 
Lo 
EoEr 
= region 2 3.3 
N= tant) (region 2) (3.3) 


the reflection coefficient [ and the transmission coefficient t are given by 


_ 12 COS 0; — 119 COS 8; (3.4) 
To COS O; + No COSA; 


272 COS 0; 


Se 3.5 
No COSO; +N COSO; as) 
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These preceding equations are known as Fresnel's formulas for parallel 
polarization [Ref. 3]. The phase constants for the material on each side of the 


boundary can be written as 


Bo = @yHoFo (3.6) 


where w =2nf (f=frequency in Hertz). The term tand in Equation (3.7) is the 
loss tangent of the material in region 2. The attenuation constant for a lossy 


dielectric material is given by 


eas a HOFOet} 4/14 (tand) 1] | (3.8) 


Refering to Figure 9, given a distance to the boundary of Rj and the distance 


from the boundary to the point of interest Ra, the electric field is given by 


tT 


— Tt _ -(O%gRo+i(BQRy+BaRo)) 
~ 4n(Rr) er) 


Pa 


E(Rr) 


The transmission coefficient t accounts for the loss due to the reflection at the 
interface, the attenuation constant a» accounts for the loss due to the material 


dielectric and loss tangent, and the phase constants B, and Bs» account for the 


changing phase with respect to a free space path. The total transmit/receive 
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distance of Ri and Ra is denoted as Rt. The value E(Ry) represents the total 
transmitted field of an oblique planar wave incident at an angle 9). 

3. Plane Waves at Oblique Incidence on a Double Interface 

The transmission and reflection of a plane wave at oblique incidence to a 
two boundary interface is shown in Figure 10. The regions on both sides of the 
lossy dielectric panel are assumed to be free space. In this case the 
transmission of the plane wave must include the multiple reflections between 
the two faces. The transmission through each boundary is governed by the 
equations previously shown for the single interface. In the two interface 


problem, refraction causes the apparent location of the source as viewed in 
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Figure 10. Plane wave at oblique incidence - double interface 





region 3 to move as shown. The original source can be replaced by an 
equivalent source which is located closer to the wall and modified to account for 
the effects of the wail. The equivalent source takes into account the 
transmission coefficients of each region based on the original source and 
geometry. The equivalent source will also take into account the phase change 
and attenuation due to the wall. The wall can be removed and the equivalent 
source will radiate in free space. In this model, multiple reflections inside of the 
wall are neglected. Note that the equivalent source method can be expanded 
to include multiple reflections. An additional equivalent source can be defined 
for each multiple reflection added. 


The transmission coefficient t, is given by 


= 272 COS O; (3.10) 
Tio COS 8; + No COS 9§; 
and the transmission coefficient to is given by 


94 = ——-e.0.wt™ 
. No COS 9; + No COS O; 


The phase and attenuation constants of the lossy dielectric are given by 


Bo =@ bata 1+ (tan 8.) +1] (3.12) 

















Qo =@ baal 1+ (tand2)? -1| (3.13) 


The electric field to the right of the second boundary can then be thought of as 


originating from an equivalent point source of the form 


(R;) = aan 97 (-aRo ti(Bo(Ry+R3)+BaRo ) (3.14) 


which radiates in free space with the wall removed. The distance to the 
boundary is Ri, the distance traveled through the material is R2, the distance 
traveled outside the second boundary is R3, and the equivalent source distance 
is Rj. The transmission coefficients account for the loss due to the reflection at 
each interface. The phase constants account for the changing phase with 


respect to the transmitting reference. 


B. PROPAGATION OF A SPHERICAL WAVEFRONT 

1. Direct Path 

The propagation of a spherical wave through a wall can be modeled by 
assuming that the radius of curvature is large and the wave front is locally 
plane. This is the standard assumption for geometrical optics [Ref. 4]. As 


shown in Figure 11, the propagation of a local portion of the wave front can be 


represented by a single ray at an angle Qj. 
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Figure 11. Locally plane sections of a spherical wave 


The maximum and minimum angles that intercept the edges of the 
receive antenna can be traced back to the point where the lines intersect as 
shown in Figure 12. This point is the location of an equivalent spherical source 
radiating in free space, and provides the geometrical optics divergence factor 
(/R;). Every ray with angle of incidence 9; will trace back through the 
intersection point. Each locally plane section corresponding to a particular 
angle of incidence is traced through the wall and the transmission coefficient 
and phase and attenuation constants calculated in each region. The rays 
coming from the equivalent source will contain the transmission coefficient and 


phase and attenuation constants from the path followed by the original rays. 
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spherical 
source 


location of apparent 
spherical source as 
viewed from region 3 





region 1 | region 2 region 3 


Figure 12. Equivalent source location 


Using this approach the equivalent source includes the loss due to 
spreading along a single straight path Rj with no wall present as shown in 
Figure 13. The radiation pattern of the equivalent source Ej is given by 
Equation (3.14). Each variable in the equation is a function of the incidence 
angle. The transmission coefficient, phase constant, attenuation constant, and 
distances Ri, Re, and Rg are based on the ray paths from the original source. 
The distance R; is the distance from the equivalent source to the point on the 
antenna grid, as indicated in Figure 13. The incident angle increment is chosen 
such that the spacing between each point on the receive aperture plane is 
approximately one tenth of a wavelength (.14). The total received field is found 


by integrating the field over the aperture [Ref 5]. 


2 1 





_ edge ray receive 
original source aperture 


plane 


equivalent source 





Figure 13. Equivalent source and antenna integration grid 


2. Geometric Relationships for the Direct Path 

Using the ray tracing concept, a building wall propagation loss model 
was formulated. In this section, formulas for the quantities appearing in 
Equation (3.14) are derived. As shown in Figure 14, the source is located at a 
height h above the floor and at a distance di from a wall of thickness L. The 
receive location de is assumed equal to d;. The regions on either side of the 


wall are free space. A typical direct ray is shown in Figure 14. 
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as 


Figure 14. Direct rays 


Given that 6; is known, the distances R;, Re, Rg and Rj can be calculated 


d 
R, =—L 3.15 
! cos 9; \ 
Ae -—_;— (3.16) 
COS@ sin (£2 sing;) 
2 
d 
Rg =~ 3.17 
7. cos 6; \ 
d,+L+ds.-di 
R, = —-—_+—+ 3.18 
i c08 6, (3.18) 
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The phase and attenuation constants of the wall are given by Equations (3.7) 
and (3.8). The transmission coefficients as a function of 8; are given by 
Equations (3.10) and (3.11). The intrinsic impedance of free space n. and the 
material no are given by Equations (3.2) and (3.3). Thus, all of the quantities in 
Equation (3.14) are known, and the electric field at the aperture can be 
computed. Using a subscript d to denote the direct path, Equation (3.14) can be 


written as 


Ey (Ri) = = (Ry+Fg)}+BaR)) 


(3.19) 
The incident angle is determined by discretizing the receive aperture into 


a grid of points. Each location within the aperture of diameter d can be 


expressed as the magnitude of Axm and Ayn, shown in Figure 15. 
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Figure 15. Antenna grid - incident angle 
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The magnitude pj is used to calculate the incident angle 0;__ 


VA 
6,4, =tan- (aetal- an int Ma (3.20) 


d,+bL+d5-d d,+b+do5-d 


The equivalent replacement source given by Equation (3.19) will create a field 
at each of the points on the antenna grid. In determining the received field for a 
vertically polarized antenna, only the x component is required. The E, field 
component of an x-directed electric dipole Io! in the far-field can be expressed 


in terms of spherical coordinates as 


_ —jknI, le! Ik 


E 
x ~ Anr 


= lie sin? @cos oc) (3.21) 


Therefore, the electric field of the equivalent source is normalized to include 


only the polarization factor (1-sin2@cos#) for an x directed infinitesimal dipole, 


(1- sin® @cos* o)tyt2_.-(aaRoti(Bo(Ry+Rg}+B2Ro)) 


4n(R) (3.22) 


Eq(Rj) = 


Equation (3.22) is the equivalent spherical source with no wall present for the 
direct path of propagation. 

3. Reflected Path 

The floor reflection is included using an equivalent image source. The 
location of the original image is directly below the source. Since the floor is not 


a perfect conductor, the image is not of equal intensity, but reduced by a factor 
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T,. With the image present, the floor is removed while the wall remains. This 
step is not entirely accurate, because it assumes that rays hitting the wall below 
floor level pass through the wall medium. However, for the geometry under 
consideration (equal antenna heights and equal distances from the wall) the 
approximation is accurate. 


The original source is traced to the receiving antenna to determine the 


maximum and minimum angles of incidence of the reflected rays, 9 and 


Fmin 
Or: rhe rays corresponding to the maximum and minimum angles of 
incidence are traced back to the point of intersection, as shown in Figure 16. 


This is the location of the equivalent floor reflected source. The geometrical 
relationships are similar to those for the direct rays. The distances Ry, Ro, and 





Rg, are 

dy 
R, = 3.23 
+ cos 0, > 
Ro, = =. re —, (3.24) 

COSY sin“ 1(E2 sino,) 
Bo 

do 
R3 = S25 
e cos0, \ 


The equivalent image source is determined by finding the point where the 


straight path lines Rr (6. ) and Rr (0, ___) intersect. Due to the non-symmetry, 


min Tmax 
the lines will intersect off of the z axis. An equation for each line is used to find 


the intersection point in dy, and d;,. The straight path length Rr for each reflected 
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Figure 16. Equivalent image source - reflected waves 
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angle 6, is given by, 


R= ———_-—_. (3.26) 


The phase and attenuation constants due to the wall are the same as for the 
directs rays given by Equations (8.7) and (3.8). The image beneath the floor will 
have less intensity than the source because the floor is not a perfect conductor. 
Thus, a reflection coefficient must be included to account for the reduced 


intensity of the floor image 


_ Tg COS O4, — No COS OF 


= (3.27) 
Te COS Ot, + No COS OF 


I, (8, ) 


where nf Is the intrinsic impedance of the floor given by 


| Ho 
= ot (3.28) 





T) = 
(1—jtands) 
and ©; and @;, are given by 
OF = (7% - 6} (3.29) 
= sin-1{ Bo/ gj 
O;, = sin (Pog, sine | (3.30) 
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The electric field resulting from the floor reflected rays which includes the loss 


due to floor reflection is 


a OY Oy sap; +)(Bo(Ry, +R, )+BaR2_ )) 


aa (3.31) 


E, (Rr) 


An equivalent image source with a pattern given by Equation (3.31) radiating in 
free space will create the same field at each of the points on the antenna grid as 
the original problem. Again, as in Equation (3.22), the source is modified to 
include only the polarization factor (1-sin¢@cos) for an x directed infinitesimal 


dipole. The total field resulting from the reflected rays is given by 


(1- sin? @cos* )I';t4, to, 9 (tae, +(Bo(Fy, +g, +B Re, ) 


an(R.) _ (3.32) 


E-(R,) 7 


4. Beam Shaping 
The effect of the antenna pattern is easily included in this model. The 


3dB beamwidth, which varies with wavelength (A) and antenna diameter (d), is 


approximately 


654 
8348 =—,— (deg) (3.33) 


for uniform excitation. To account for the decreasing beamwidth with increasing 


frequency, a cosine taper function was used which is raised to the nth power 


XN 
cos"(6 = cos" (65 ———_) =. 707 34 
(838) ( a 1a0. (3.34) 
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where the exponent is given by 


_ __log(.707) (3.35) 


Thus, the direct and reflected fields are Eycos"6; and E, cos" 6, where Ey and 
E, are defined in Equation (3.22) and (3.32), respectively. The direct and floor 
reflected fields are now a function of antenna beamwidth, which varies with 
frequency. The beamshape has very little effect on the direct term, but 
significant effect on the floor reflected term. 

5. Total Received Field 


The complex field of the direct and reflected rays are summed to 


‘determine the total resultant field Ey. From aperture theory [Ref. 6] the total far 


field received is an integral of the aperture fields 


. ee IKl ras 
alee —- J) IEA hid ods, S28) 
SA 


where Sa is the aperture surface, k is the incident ray propagation vector and 
d the position vector to the point of integration. Neglecting constants and 
writing Equation (2.36) in terms of previously defined quantities for a discretized 


grid 


Et = > ¥(Edan +E tin |? AxAY (3.37) 
m nA 
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C. MODELING RESULTS 


1. Overview 

The spherical wave propagation model was programmed in MATLAB (a 
code listing is provided in Appendix B). If a single frequency is input, the 
program calculates and plots the magnitude and phase of the direct, reflected, 
and total fields on the surface of the antenna and provides a loss in GB to that of 
free space for the user specified wall parameters. If more than one frequency 
point is input, the program calculates the loss in dB for each frequency and piots 
the resulting loss versus frequency over the band specified. The program 
variables specified by the user are: 


di, d2 = distance from antenna to wall (ft) 


L = thickness of wall (ft) 

h = height of antennas from floor (ft) 
e2 = dielectric constant of wail 

It = loss tangent of wall 

ef = dielectric constant of floor 

tf = loss tangent of floor 


2. Comparison of Computed and Measured Data 

Using the wall loss data from the measurements described in Chapter lH, 
Section B, comparisons were made to the data computed using the model in 
Chapter III, Section B. The wall characteristics are summarized in Table 2. The 
electrical characteristics of the walls were not known precisely, so parametric 
calculations were performed. The quantities e«, and tand were varied until the 


best match with measured data was obtained. Comparisons of the measured 


31 








and computed insertion loss data are given in Figures 17 through 21. A 


summary of the dielectric and loss tangents are given in Table 3. 
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Figure 17. Insertion loss of 25" glass pane 
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Figure 18. Insertion loss of 1.75" wood doors 
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Figure 19. Insertion loss of 5.5" wall w/ chalkboard 
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Figure 20. Insertion loss of 10.75" cement wall 
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Figure 21. Insertion loss of 1.75" metal doors 
The outputs for the single frequency calculations are provided in 


Appendix C. These plots include the calculations of field values over the grid of 


points for the direct, reflected, and combined terms. 
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Table 3. Summary of propagation loss measurements and calculations 



















ue 


























3. Summary 
A mathematical model was developed to approximate the propagation of 
a spherical electromagnetic wave front through a wall. Discrepancies between 


the measured and computed data can be attributed to the limitations of tne 
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model and possible inaccuracies in the measured data. The assumptions 


inherent in the model are as follows: 


1) the polarization of the source and receive antennas is vertical, 

2) surface roughness of the wall is neglected and the faces are parallel, 
3) loss tangent and dielectric are constant throughout the wall, 

4) a frequency independent dielectric constant and loss tangent are 
used in the equations for impedance, phase, and attenuation, and 


5) multiple reflections within the wall are neglected. 


Based on the data in Figures 17 through 21 several conclusions can be made. 
First, it does not appear that the effects of mulitpath are significant. The ripple at 
lower frequencies is at most 4 dB peak to peak for low loss materials such as 
wood (Figure 19). For lossy materials the multipath effect is eyen less 
noticeable. Secondly, one can conclude that this simple model provides a 
relatively accurate estimate of wall losses if the dielectric constant and loss 
tangent can be determined from visual inspection of the wall. 

With regard to the MRPV, the smaller wall losses at lower frequencies 
indicate that the 2.45 GHz scientific and commercial frequency band is 
preferable to the 5.5 GHz band. However, the antenna design also enters into 
the frequency selection tradeoff. As the frequency decreases, the gain 
decreases for a fixed antenna size. This effect favors the higher frequency. The 
wall loss is probably the more important of the two, because the beam may pass 


through several walls. 
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iV. ANTENNA DESIGN 


A. BACKGROUND 

The transmission of beamed microwave power for the MRPV is designed 
such that a substantial amount of pitch and roll can be tolerated without 
significant power collection loss. This is accomplished with a highly directional 
transmitter and a near omni-directional receiver. Therefore, the reception of 
power for the MRPV is based on an antenna design that provides maximum 
coverage in space. This coverage is attained through a design that 
approximates the far-field pattern of an ideal dipole. 

1. Configurations 

In determining the antenna design for the MRPV, several configurations 
exist for both horizontal and vertical polarizations. Given consideration of the 
effect of ground returns on the total received field, the polarization of the 
antenna is of importance. The multi-path effect that is present when operating 
near the ground can cause an undesired cancellation of the direct and reflected 
paths of the transmitted microwave beam. This would result in a reduction in 
the total received field at the antenna. A wave with horizontal polarization has a 
large ground reflection coefficient while a wave with vertical polarization has a 
smaller ground reflection. Thus multi-path favors vertical polarization. 
However, multi-path is only of secondary concern because the transmit antenna 
beam is very narrow. The major reason for using vertical polarization is the 
most efficient and compact antennas that conform to the vehicle shape are 
vertically polarized. Two such antenna configurations for the MRPV are shown 


in Figure 22. 
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Figure 22. Possible antenna configurations 


2. Circumferential Slot | 

The preferred configuration for vertical polarization is the circumferential 
slot antenna. As shown in Figure 22, the design uses direct excitation of the 
body with feed points on each of the insulated body halves. The primary 
advantage of the slot antenna over the array of vertical dipoles is the simplicity 
of the feed. The vertical array of dipoles would consist of at least four dipoles 
requiring some feed circuitry. The complexity of the design would also 
increase based on the number of dipoles chosen in the array. The radiation 
pattern of the dipole array will also be affected by the mutual coupling between 
the dipoles themselves as well as with the vehicle body. The results of the 
propagation loss study (Chapter Ill, Section C) indicate that the 2.45 GHz 
scientific and commercial frequency band is preferable to the 5.5 GHz band. At 
2.45 GHz the circumferential slot with a cylinder length of approximately 6 cm 
corresponds to an excitation of (1/2)A. At 5.5 GHz the slot antenna with a 


cylinder length of 6 cm corresponds to an excitaion of (1)2. 
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Although the electric field pattern is optimized for maximum coverage (i.e. 


has a dipole pattern) at frequencies where the body is less than a wavelength, 
the efficiency of the antenna favors the lower frequency. Antennas that are 
electrically small have a large input reactance. In the case of the MRPV, this 
problem is avoided because, although the antenna is electrically short, it is not 
electrically thin. The antenna impedance is important in the transfer of power 
from the receiver to the battery. Since the final dimensions of the MRPV are not 
yet chosen, the design of the circumferential slot is centered around a full 
wavelength excitation. Due to the limitations of testing in the NPS anechoic 
chamber, a scaled version of the antenna prototype was designed to operate at 


X-band, whereas the actual antenna will operate at S-band. 


B. SIMULATION 

1. Method of Moments 

The calculations of far-field patterns for the prototype antenna designs 
were done using the PATCH computer code developed at the Sandia National 
Laboratories (Ref. 7]. PATCH is a frequency domain electromagnetic scattering 
code based on a method-of-moments solution to the Electric Field Integral 
Equation (EFIE). The objects in PATCH are modeled by planar triangular 
patches, which easily conform to surfaces and boundaries of general shape 
and allow variable patch densities over the surface of the object. Although the 
EFIE has been formulated in terms of a scattering problem, the PATCH code 
has the ability to treat both plane wave and voltage source excitations. 
Therefore, PATCH is capable of handling antenna problems if the terminals 
across which a voltage drop can be defined are identified. The voltage drop 


can be represented by cutting a thin slot along the edge of two adjoining 
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patches to create a potential difference between the two triangular patches. 
This is shown in Figure 23 where the two adjoining triangular patches, A and B, 
are specified to have a potential difference of Vg with a surface current Je 
across the slot. To model a structure effectively, the density of the triangular 
patches must be sufficient to accurately represent variations in surface 
geometry and surface current density. The edge lengths of each triangular 


oatch used are 1/10 of a wavelength. 


- Vg + 


IS INAENTAAAN 


antenna geometry surface current orientation 


NES —-09) TFANTAI 


potential difference at antenna terminals 





Figure 23. Voltage source excitation with PATCH 


2. Antenna Model 

Two models are constructed for the circumferential slot antenna: 1) an 
edge fed slot antenna, and 2) a center fed slot antenna. Both models were 
used to determine the effects of feed point location on the field pattern. in the 
first case, a triangular mesh structure is built for two unattached cylinders each 
closed on one end. The edge fed antenna model is shown in Figure 24, where 
the upper and lower cylinders are connected on the edge by a quadrilateral tab. 


The tab, as shown in Figure 25, serves as the feed point for the antenna 
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terminals. The units have arbitrarily been chosen so that the cylinder length is 
500 m, and therefore a frequency of 0.6 MHz was used in the computation. 


This choice results in an overall antenna length that matches the desired (1)A. 


MAK DIN 
NNN 


A 
Sa 





Figure 24. Edge fed slot antenna - isotropic view 
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Figure 25. Edge fed slot antenna - profile view 


The center-fed antenna model is shown in Figure 26. In this case, the 
upper and lower cylinders are closed on each end. A quadrilateral tab 
structure, which serves as the feed point for the antenna terminals, is located on 
the center point of each closed end as shown in Figure 27. The tab structure is 
constructed with a symmetrical center feed to provide a current with matching 
phases along the perimeter of the cylinders. Again, the excitation frequency is 
scaled relative to the overall antenna length to obtain the desired size in 


wavelengths. 
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Figure 27. Center fed slot antenna - profile view 











3. Results 

The normalized computed patterns for the edge fed siot antenna are 
provided in Figures 28 and 29 for the E and H plane far field patterns, 
respectively. If the vehicle is flying straight and level, the E-plane (or vertical 
plane) is perpendicular to the ground; the H-plane (or horizontal plane) is 
parallel to the ground. As shown, the location of feed points results in a non- 
symmetric field pattern with nulls in both planes. This is caused by the surface 
currents being out of phase along the perimeter of the antenna. Therefore, the 
location of the feed point cannot be on the rim of the body, but must be centered 
on each cylinder interior cap to provide equal phasing of the surface currents. 
The normalized computed results for the center fed slot antenna are provided in 
Figures 30 and 31 for the E and H plane far field patterns, respectively. With a 
center feed point, the surface currents are flowing radially outward from the 
center of the body and arrive in equal phase along the perimeter of the body. 
As shown, this results in a symmetric field pattern in the E plane, and a near 
constant field pattern in the H plane. As a result of the PATCH simulation, a 


center fed circumferential slot antenna is chosen for the prototype design. 
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Figure 28. E-plane of edge fed slot antenna 
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Figure 29. H-plane of edge fed slot antenna 
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Figure 30. E-plane of center fed slot antenna 
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Figure 31. H-plane of center fed slot antenna 
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C. ANTENNA PROTOTYPE 

A prototype of the center fed circumferential slot antenna was designed 
and tested. The antenna, depicted in Figure 32, consists of two cylindrical 
copper (Cu) caps attached to a 50 ohm RG 141 semi-rigid coaxial cable. The 
center and outer conductors are soldered to the caps, which are then closed on 


each of the ends with a Cu film to enclose each cylinder. Since the antenna is 





Figure 32. Prototype antenna 


symmetrical in nature, the currents excited from an incident wave will also be 
symmetrical (or balanced). However, when an incident wave Is launched down 
a coaxial line to a symmetrical antenna, a current will flow back on the outside 
of the outer conductor which unbalances the currents on the antenna and 
transmission line. The unbalanced currents will degrade the ability of the 
antenna to radiate and increase its input VSWR. In order to balance the 
currents on the coaxial line, a quarter wave balun is connected to the antenna. 


The balun, which is a contraction for "balanced to unbalanced," suppresses the 
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outside current which is causing the imbalance. The balun transforms the 
balanced input impedance of the slot antenna to the unbalanced coaxial line 
such that there is no net current on the outer conductor of the coax [Ref. 5]. A 


picture of the prototype slot antenna is shown in Figure 33. 





Figure 33. Circumferential slot antenna - prototype 
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V. ANTENNA MEASUREMENTS 


A. ANTENNA PARAMETERS 

The antenna parameters were measured using an HP 8510C network 
analyzer. The reflection parameters of the network were measured to 
determine the impedance, return loss, and Voltage Standing Wave Ratio 
(VSWR) of the slot antenna. Using the HP 8510C, a scattering parameter 
- measurement was made to determine the S parameters S;; and Syo of the slot 


antenna. The measurement set-up is depicted in Figure 34 [Ref. 8]. As shown, 


HP 8510C incident transmitted S21 
a 


Do 


reflected s,, 





Figure 34. Scattering parameter measurement 


the S11 parameter represents the reflection coefficient [, and the Sy 
parameter is the forward transmission coefficient t. The measurements were 
taken from 9.5 to 10.5 GHz. The return loss is extracted from the S parameters 
and provided as a plot in Appendix D. A smith chart of the antenna impedance 


is included as Figure 35. The complex antenna impedance varies with 


frequency from Zosauz = 23.2-§9.3 to Zigs5GHz = 98.6+j36.8. The optimum 
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impedance for maximum transfer of power with a 50 ohm coaxial line is 
Z,=50+j0. At a frequency of approximately 10.3 GHz, the complex 
impedance consists of a real 50 ohms with an inductive reactance. Thus, by 
tuning the antenna in frequency and matching it with a conjugate reactance, the 
antenna can be matched to the coaxial transmission line. The measured return 
loss at the full wavelength excitation of 9.9 GHz is -9.7 dB. Therefore, the power 


reflection coefficient is equal to 


P|? = 10-089/10). — 49-9719) = 0.107 | (5.1) 


Thus about 89% of the power delivered to the antenna is radiated, and 10% is 
reflected. This is a mismatch loss of about .45 dB. The VSWR of the antenna at 


9.9 GHz is 


1+|0|_ 1+|0.327| _, 9, 


pel se = 1. 5.2 
{|r| ~ 1-|0.3271 me 


VSWR = 


B. PATTERN MEASUREMENT 

The far-field patterns of the prototype slot antenna were measured in the 
NPS anechoic chamber. A general layout of the X-band anechoic facility is 
provided as Figure 36. As indicated in Chapter IV, Section A, the antenna 


prototype is designed at a scaled version to operate at X-band. 
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Figure 36. NPS X-band anechoic chamber 





1. Test Set-Up 

The test set-up for the far-field pattern measurements is shown in Figure 
37. The HP X890 transmitting horn is oriented for horizontal transmission for 
the electromagnetic wave. For an E plane measurement, the slot antenna is 
centered over the pedestal with a horizontal orientation. The transmitting horn 
is oriented with a horizontal E plane and the pedestal then rotated with 
measurements taken for the specified angular increment. The H plane 
measurement is taken in the same manner, with the transmitting H plane 
horizontal, and the slot antenna centered over the pedestal with a vertical 
orientation. A 386 PC with the test equipment software Lab View is used to 


operate the test equipment. 
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Figure 37. Far-field pattern measurements 


2. Far-Field Measurements 

The normalized far-field patterns for the E and H planes of the slot 
antenna were measured and are provided in Figures 38 and 39, respectively. 
The E plane measurement from 0° to 180° is shown to have a half power 
beamwidth HPBW=102°. The portion of the pattern from 180° to 360° is 
affected by the exposure of the balun and cable blockage of the incident field, 
thereby corrupting the pattern. The H plane measurement contains a ripple with 
a maximum value of 8 dB. Figure 38 indicates that the pattern is also affected 
by the exposure of the balun and shadowing of the incident field by the feed 


cable. 
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Figure 38. Measured E 
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Figure 39. Measured H-plane pattern 
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C. RF TO DC POWER CONVERSION 

The final section of measurements deals with the conversion of the 
microwave signal to usable dc power. The device that accomplishes this is 
known as the rectenna and consists of the slot antenna and a diode rectifier. 
The antenna collects the incident field and delivers microwave power to the 
feed point. The rectifier converts the microwave power to a dc signal. In order 
to determine the rectenna's ability to convert power, the efficiencies of the 
antenna and rectifier are estimated. 

1. Antenna Efficiency 

The efficiency of the slot antenna is determined by measuring the gain of 
the antenna and comparing it to the directivity. The antenna efficiency is 


defined by 


G 
=— o) 
: D pe 


where G is the antenna gain and D is the directivity. The directivity represents 
the power gain an antenna would have if all received power appeared as 
output power. The gain of the slot antenna is measured using a constant power 
transmitter and a reference antenna with known gain. Using the known 
standard gain of the reference antenna, the absolute maximum power gain of 
the slot antenna is calculated from the differences in the measured receive 
powers. The ratio P;/P, is the increase in received power from the test antenna 
over that of the standard gain antenna. Therefore, the measured gain of the test 


antenna is given by [Ref. 5] 
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Pp 
Gy =5-Gs : (5.4) 


S 


The Microline 56x1 antenna is used as the standard gain horn antenna with a 
specified nominal gain Gs = 17 dB at X-band. The measured power gain of 
each antenna for a constant transmitted power is Pt = -38.25 dB and Ps = -21.75 


dB. Writing equation (5.4) in terms of dB, yields 


G, (dB) = P;(dB)-P,(dB)+G,(dB) (5.5) 


The gain of the slot antenna is determined from 
Gt = -38.25 dB - (-21.75 dB) +17dB=.5dB . (5.6) 


This measurement of gain was taken at @ = 90° in the E plane and o =0° in the 
H plane. At o=0°, a gain of -4.7 dB is measured in the H plane with an 
average gain of -3.7 dB. The accuracy of the measurement is = .5 dB. 


Therefore, the gain of the slot antenna is now given by 
Gt = .5dB+(-3.7-(-4.7))dB=.5dB=154.5dB . (5.7) 


The directivity of an antenna is the maximum value of directive gain given by 


Um 


Uave 


D= (5.8) 
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where Uaye is the average radiation intensity defined as 


P 
Uave = ze |JU(6,9)d0 = = (6.9) 





and U,, is the maximum radiation intensity given by 


Un = — > (5.10) 


where the term P; is the total radiated power and Q, is the antenna beam solid 


angle. The antenna beam solid angle is the solid angle through which all the 
power would be radiated if the power per unit solid angle (radiation intensity) 
equaled the maximum value over the beam area. The directivity can now be 


defined as 








= =— | (5.11) 


The antenna beam solid angle is given by 


Qn = ffIF(@.o)P aa (5.12) 





where the term IF(@, 0) is the normalized field pattern replaced by two fourth 
order polynomials f,;(6) and f,o(8@) which approximate the E plane data of 
Figure 38 in a least squares sense. The approximation is shown in Figure 40 
where the normalized E plane data is linearized and taken from 0 to 2x. The 


polynomial f,;(6) exists from 0 to m given by 


fn1(8) = -.00640% —.02626°% -.09996* +.73680—.0068 (5.13) 
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theta (rad) 


Figure 40. Polynomial approximation to E plane data 
and the polynomial f,2(8) exists from m to 2x given by 
f.9(@) = -.1998(y)* + 1.3188(y)° —2.7836(y)* + 1.977 Uy)-.1237 (5.14) 


where the term y=9-2. The beam solid angle is now written as 





fax max 


Tt 1 21% 
o, = { f'n() sinededo =( ( 121! sin eaedo + fna(9-) sin gdedo (5.15) 
A 
f 
00 9.0.74 rt 0 


where each polynomial is normalized by the maximum value fmax. With work, 


the integral of Equation (5.15) results in 


Qa = 04+ Q2 = 4.924+1.08 = 6.0 Sr . (5.16) 
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Using the calculated beam solid angle, the directivity of the slot antenna is 


given by 


D= 2% = 2.09=3.21 dB . (5.17) 


The efficiency of the slot antenna is now calculated with the measured gain and 


directivity as 


G, — 10(18/10) 


> “S2TT0) =,675 =67.5% - (5.18) 


e= 


However, there is a = .5 dB gain measurement error. Therefore the efficiency 
could be as low as 60% or as high as 76%. 

2. Rectenna Efficiency 

The combined capture and rectification efficiency of the rectenna can be 
estimated using the antenna efficiency and a known rectifier efficiency. The 
efficiency of diode rectifiers have been measured as high as 80% [Ref. 9]. 
Assuming a rectifier efficiency e-=.80 as demonstrated in earlier experiments, 
the maximum overall power efficiency for the circumferential slot used as a 


rectenna device Is 


Cowr = Eslot sd i — (0. 76)(. 80) = O. 608 = 61 % : (5. 1 9) 
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Vi. CONCLUSIONS 


A. SUMMARY 

A propagation loss study and antenna design were presented for the 
MRPV. The concept of beamed microwave power transmission was shown to 
be limited by conversion efficiencies and propagation loss due to spreading of 
the electromagnetic wave and any blockage encountered. The reception and 
conversion from beamed microwave power to de is performed by a device 
known as the rectenna which consists of an antenna and a rectifier. 

The insertion loss of typical building walls was measured using the 
system described in Chapter Il, Section A. The walls tested consisted of low 
loss materials such as wood doors to high loss materials such as metal covered 
doors. An approximate model was developed for the propagation of a spherical 
electromagnetic wave front through a building wall of known material 
properties. The model is based on the method of images and geometrical 
optics. Using locally plane wave approximations and ray tracing techniques, 
the signal strength of each plane wave was determined at the receive antenna. 
Combining the direct and reflected paths, the total resultant field strength was 
calculated through summation and complex integration. The model was shown 
to provide relatively accurate results. The measurements and modeling were 
performed over a band from 2 to 6 GHz. Two primary frequency bands of 
interest were looked at for overall system design: 1) the 2.45 GHz scientific and 
commercial band, and 2) the 5.5 GHz satellite uplink frequency band. The 
smaller wall losses at the lower frequencies indicate that the 2.45 GHz 


frequency band is preferable to the 5.5 GHz band. The gain of an antenna with 
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a fixed size decreases with lower frequency, thereby favoring the higher 
frequency. However, the wall loss was determined to be the more important of 
the two because the beam may pass through several walls. Thus the 2.45 GHz 
band (S-band) was used for overall system design. 

The transmission of beamed microwave power for the MRPV requires a 
highly directional transmitter and a near omni-directional receiver. The design 
of the receive antenna provides maximum coverage in space. Two possible 
configurations were presented for a vertically polarized antenna - the 
circumferential slot and the vertical array. Due to the simplicity of a single feed 
point, the slot antenna was the selected approach. Using the PATCH computer 
code, simulations were performed on two slot antennas - the edge fed slot, and 
the center fed slot. The far-field patterns were calculated and the center fed slot 
antenna was shown to be the preferred design. A prototype of the 
circumferential slot antenna was designed, built and tested. Antenna 
parameters and far-field patterns were provided, with an estimated antenna 
efficiency between 60 and 76%. The sources of loss are feed point mismatch, 
feed cable loss, and ohmic losses due to the imperfectly conducting antenna 


surfaces. 


B. FUTURE CONSIDERATIONS 


The continuing research efforts for the MRPV should concentrate on the 
areas of microwave to dc power conversion (rectenna), improvements to 
propagation loss model, and transmitter design. The microwave to dc power 
conversion includes the design of a diode rectifier with sufficient conversion 
efficiency to supply power for operation of a dc motor. The diode rectifiers 


looked at during the course of this study were the HP 8471E series. These 
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diodes were designed as a measuring device and were not intended for power 
conversion. Improvements to the propagation loss model can be made in the 
areas listed as assumptions in Chapter Ill, Section 3. The inclusion of mulitple 
reflections and frequency dependent dielectric constants and loss tangents will 
improve the accuracy of the results. Finally, with the results of the propagation 
loss study and rectenna design, the performance requirements of a highly 


directional transmitter and antenna can be specified. 
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APPENDIX A 
WALL LOSS DATA 
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Figure 41. Loss measurement of .25” glass pane, 2 to 2.9 GHz 
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Figure 42. Loss measurement of .25” glass pane, 2.9 to 6 GHz 
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Figure 43. Loss measurement of 1.75” wood doors, 2 to 2.9 GHz 
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Figure 44. Loss measurement of 1.75” wood doors, 2.9 to 6 GHz 
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Figure 45. Loss measurement of 5.5” wall w/ chalkboard, 2 to 2.9 GHz 
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Figure 46. Loss measurement of 5.5” wail w/ chalkboard, 2.9 to 6 GHz 
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Figure 47. Loss measurement of 10.75” cement wall, 2 to 2.9 GHz 
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Figure 48. Loss measurement of 10.75” cement wall, 2.9 to 6 GHz 
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Figure 49. Loss measurement of 1.75” metal doors, ¢ to 2.9 GHz 
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Figure 50. Loss measurement of 1.75” metal doors, 2.9 to 6 GHz 
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APPENDIX B 
MATLAB PROGRAM 


% Thad B Gibson 
%o 
% RAY4.m - BUILDING WALL PROPAGATION LOSS MODEL 

% 

% This program models the propagation of a spherical electromagnetic 

% wave front through a wall. The program will calculate an equivalent 

% replacement source for the DIRECT and an equivalent image of the 

° REFLECTED portions of a spherical wave. If a single frequency is chosen, it 
°% will calculate and then plot the magnitude and phase of the DIRECT, 

% REFLECTED, and COMBINED rays on the surface of the antenna and 

% provide a Loss in dB to that of free space. If morethan one frequency point Is 
°%% chosen, it will calculate the Loss in dB for eachfrequency and then plot the 
% resulting Loss vs. frequency. 

% 

% 

clear; clg; format compact; 

% 

% 


% rrr rerrrrrrrrrrrrrrrrrr rrr rrr Trt ttt hte he od hed kd kk iaksialieialaieiaheieieheieelel 
o 


% ** INPUT - SET UP ALL INPUT VALUES si 
% Fee FIT TOI IIT TT TI IOI IORI TF TOTO TA TOIT IT TFT TE IT IIT TAS TAKIN SAK EIA IIIS ASIAN IRIE 
% 
oflag=input(‘Do you want hardcopy?(y/n) ',’s’); 
fstart=input(‘Enter the start frequency (GHz) ’); 
fstop=input('Enter the stop frequency (GHz) }); 
count=input(‘Enter # of points to calculate between fstart and fstop '); 
if fstart == fstop, 


count = 1; 
otest=input('Do you want plots of phase as well as magnitude?’','s’); 
end; 


d1=input('Enter the distance from antenna to wall (ft) ‘); d2=d1; 
l=input(‘Enter length of wall (ft) ‘); 

h=input('Enter height of antennas from floor (ft) ’); 
erx=input(‘Enter the dielectric constant of wall ’); 
tandx=input('Enter the loss tangent of the wall’); 
erfx=input(‘Enter the dielectric constant of the floor ’); 
tandfx=input('Enter the loss tangent of the floor’); 
fstart=fstart* 1e9; 

fstop=fstop*1e9; 

% 

% ** END INPUT a 
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%o 
% KKK KK KKK HK HH KK KKK IKK KKK KKK KKK IKKE KEE KER KKK KEKE KEKE EERE KKKKKEEKEKKKEKKKKKKKEK 
%o *™ OUTER LOOP - SETS FREQUENCY BEING CALCULATED ~ 
% KK KKK KKH HEE KKK KEKE KRKKEREKEKKKKKKKKKKKRRKKKRKKKKKRKKKKKKKKKKKKKKKKKKK KKK KKKKK KKK KKKEK 
%o 
Yo 


frange=linspace(fstart,fstop,count); 
for freq=1:count, 
f=frange(freq); 
c=3e8; d=1.5:; 
er=erx; erf=erfx; tand=tandx; tandf=tandfx; 
incr=round((d/3.281)/(.1*c/f)); 


% 
% 
% TC REESRALEESERAEEESRSESEERASEEREARASAARERBAEEERASERERRARASERELELASZALESESAE SS ASLE LAR ERE LAS ESE LE SS 
% ** INNER LOOP - USED TO DETERMINE A LOSS FOR A GIVEN si 
%** FREQUENCY BY CALCULATING WITH WALL AND SUBTRACTING ** 
% ** FROM WITHOUT WALL (FREE SPACE) = 
% KKK KAKI KEKE K KE KEKE KEKE EKER KERKREKKKCEEKKKEKKEKRKKKEKKKKKKKKKKKKKKEKKE 
% 
% 
for loop=1:2, 

if loop == 2, 

er=1; erf=1; tand=0; tandf=0; 

end 
% 
% Calculate material constants 
% 
e0=8.84e-12; uo=12.56e-7; 
w=2"pi't; 


alphad=w*sart((uo*eo*er/2)*(saqrt(1+tand’2)-1)); 
alphaf=w*sart((uo*eo*er/2)*(saqrt(1+tandf2)-1)): 
bo=w*sart(eo*uo); 
b2=w*sart((uo*eo*er/2)*(sqrt(1+tand42)+1)); 
bf=w*sart((uo*eo*er/2)*(saqrt(1+tandf42)+1)); 
no=saqrt(uo/eo); n2=saqrt(uo/(er*eo))/saqrt(1-j*tand); 
nf=sqrt(uo/(erf*eo))/sqrt(1-j*tandf); 

% 

% Initialize arrays 

% 

Ei=zeros(incr); Er=zeros(incr); Et=zeros(incr); 
theta=zeros(size(1:incr)): 
ohi=zeros(size(1:incr)): 
thetar=zeros(size(1:incr)); 
ohir=zeros(size(1:incr)); 
dxi=linspace(-d/2,d/2,incr); dyi=dxi; 
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dxr=linspace(2*h-d/2,2*h+d/2,incr); dyr=dxi; 
% 
% 
% weer Tr rrrrr rrr rrrrrrerr rer rrr rrr PTI PEEL EE LST EPPS STE Pho hk elke 
o%** DIRECT RAYS - FIND MIN AND MAX INCIDENT ANGLES TOFIND ™ 
o% ** THE EQUIVALENT REPLACEMENT SOURCE LOCATION - diz ie 
% arn Sr n rrr rrrrrrrrrrrre rere rt rrr Err EEE SEL ELS SL LS tL eh the ee Rk ii ieieieialel 
% 
% 
acc=200; 
ang=linspace(0,pi/4,acc); 
for k=1:acc, 

y1=(d1+d2)*tan(ang(k))+l*tan(asin(bo/b2*sin(ang(k)))); 

if y1 > d/2, 

break; 

end 
end 
angimax=ang(k-1); 
angimin=-angimax; 
% 


o ] 
iz=d1+l+d2-((d/2)/tan(angimax)); 

% 
% 
% ** END DIRECT RAY 2 
Of, mnncwenccn nnn nnn nnn nae enen nese ene c ena nneren nnn wennecnse cnn nnnccsennnanennnanascrenancaccsencnccnsessare 
% 
% 
% 
% 
% SSS rrr rrrrrrrrrrrrrrrrrrcrroL ECL ELE LILES ELSA SPCR eR Le a ike lll 
%** REFLECTED RAYS - FIND MIN AND MAX INCIDENT ANGLES TO a 
%** FIND THE EQUIVALENT IMAGE SOURCE LOCATION - drx,drz ae 
% rere rrr rrr rrr rrr rrrrereorrcrrrcrre errr ti eT ESS EST FP Lec LF eC tPF hk ee Rlel 
% 
% 
acc=400; 
ang=linspace(0, pi,acc); 
for k=1:acc, 

y2=(d1+d2)*tan(ang(k))+I*tan(asin(bo/b2*sin(ang(k)))); 

if y2 > 2*h-d/2, 

break; 

end 
end 
angrmin=ang(k-1); 
%o 
for k=1:acc, 


y2=(d1+d2)*tan(ang(k))+I*tan(asin(bo/b2*sin(ang(k)))): 
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if y2 > 2*h+d/2, 


break; 
end 
end 
angrmax=ang(k- 1); 
% 


% Calculate drmin and drmax -the min & max deviation (z) from the original 
% source to the image 

% 

drmin=d1+d2+1-((2*h-d/2)/tan(angrmin)); 
drmax=d1+d2+I-((2*h+d/2)/tan(angrmax)); 

% 

% The relfected rays do not trace back to the z axis, so must determine 
% where they cross in (z,x) coordinates. Calculate equation of line for 
% max and min angles and find intersection in (drx,drz) 

% 

bmin=-tan(angrmin)*drmin; 

bmax=-tan(angrmax)*drmax; 

mmin=tan(angrmin); 

mmax=tan(angrmax); 

% 

% Now find the point where the two lines cross (drx,drz) 
drz=(bmin-bmax)/(mmax-mmin); 

drx=mmin*drz+bmin; 

% 

% Now find distance from point (drx,drz) to center of antenna 
Re=sart((2*h-drx)42+(d2+l+d1-drz)*2); 

% 


% 
% ** END REFLECTED RAY aia 
Of, nawennnannnncecnennennnacneececnnnneneensccnaeseccnnenannensneenenenneneneeaeneneeneersananaenaecnerne= 


% renner rrr Terr rrrrrrerroeecrCeceCCLCCC CC SER EE LES ES ERLE RSE ARES EL CSAS ee 


% ** EFFECT OF A TYPICAL ANTENNA BEAMWIDTH WHICH VARIES . 


%** WITH FREQUENCY ‘a 
% rer rrrrrrrerr orrrrrercrecroereecrrcrccCSSOOS££S ESSE REL ESE RE RELEASE RAR RASS SD ee Rel 
Yo 
% 


% Now account for beamwidth of antenna which is approximated by (65*(c/f)/d) 
p=log10(.707)/log 10((65*(c/f)/d)*(pi/180)); 
% 


Ye" END OF EFFECT OF BEAMWIDTH i 
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rrr Sr rrr rr rrrrrrrrrrr eer r errr TST STL TPP Phe Lk eh he Reiki keke keeleieie Rel Relelel 


% 
%** DIRECT RAYS - FIND Ei FOR EACH GIVEN ANGI, THEN STORE ois 
%** EACH VALUE OF Ei INA SQUARE MATRIX = 
% rrr rrr rrrrrrr rere rr ero rrr Err LETS PTE ST Leh he ee hk eee oki Reiki R Rei ieieieielel 
% 
% 
for m=1:1NCr, 

for n=1:iNCcr, 

rho=sart(dxi(n)42+dyi(m)*2); 
if rho < d/2, 


theta=atan(dxi(n)/(d2+l+d1-diz)); 

phi=atan(dyi(m)/(d2+l+d 1-diz)); 

angi=atan(rho/(d2+l+d1-diz)); 

R1=d1/cos(angi); 

R2=!/cos(asin(bo/b2"*sin(angi))); 

R3=d2/cos(angi); 

Ri=(d2+1+d 1-diz)/cos(angi); 

angt=asin(bo/b2"sin(angi)); 

t1=(2*n2*cos(angt))/(n2*cos(angt)+no*cos(angi)); 

t2=(2*no*cos(angi))/(no*cos(angi)+n2*cos(angt)); 

E=[t1*t2/(4*pi*Ri)]*exp(-alphad*R2-j*(bo*(R1+R3)+b2*R2)); 

Ei(incr+1-n,m)=cos(angi)4p*(1-sin(theta)*2*cos(phi)*2)*E; 
end 


end 
end 
% 
% 
% ** END DIRECT RAY CALCULATION OF Ei mn 
% ee eee nh tt 
% 
% 
% 
% 
% errr re rrr rrr rrr rrrrrrrrrrrcrrrrororerrer errr ri ett PSPS ELEC SEP EPL APA ee 
%** REFLECTED RAYS - FIND Er FOR EACH ANGLE ANGR, THEN ‘is 
%** STORE EACH VALUE OF Er IN A SQUARE MATRIX 
% ee rrr rr rr rrr rrrrrrrrreorrerrerrrreerercrroesrrreeeeLeLe ELE si cise ci hss Li cesar sks 
% 
for m=1: incr, 

for n=1:incr, 

grid=sqrt(dxi(n)42+dyi(m)*2); 
if grid < d/2, 


thetar=atan((dxr(n)-drx)/(d2+l+d1-drz)); 
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ohir=atan(dyr(m)/Re); 

rhor=sart((dxr(n)-drx)42+dyr(m)%2); 

angr=atan(rhor/(d2+l+d1-drz)); 

Rir=d1/cos(angpr): 

R2r=I/cos(asin(bo/b2"*sin(angr))); 

R3r=d2/cos(angr); 

Rir=rhor/sin(angr); 

angtr=asin((bo/b2)*sin(angr)); 

t1r=(2*n2*cos(angtr))/(n2*cos(angtr)+no*cos(angr)); 

t2r=(2*no*cos(angr))/(no*cos(angr)+n2*cos(angtr)); 

angf=pi/2-angr; 

angtf=asin((bo/bf)*sin(angf)); 

Gr=(nf*cos(angtf)-no*cos(angf))/(nf*cos(angtf)+no*cos(angf)); 

E=[Grt1 r*t2r/(4*pi*Rir)]*exp(-alphad* R2r-j*(bo*(R1 r+R3r)+b2*R2r)); 

Er(incr+1-n,m)=cos(angr)4p*(1 -sin(thetar)42*cos(phir)42)*E; 
end 


% Perr Tarr rrrrrrrrrrrrrerrrr rrr ts SPP TTT Pitt et hh eh oe eh he ake ike eeiekeieieeletel 
° 


% ** DETERMINE MAGNITUDE AND PHASE OF DIRECT, REFLECTED ™ 
%** AND TOTAL ae 
% errr rr re rrrrrrrrrrre rrr rr rest rrr eT ST tt teehee eek kiki ieiekieiteiieieielel 
% 
% | 
% Now determine the magnitude and phase across the antenna for both the 
% DIRECT and REFLECTED rays 
% 
magi=abs(Ei); phas=phase(Ei).*(180/pi); 
magr=abs(Er); phasr=phase(Er).*(180/pi); 
% 
% Determine magnitude and phase of COMBINED direct and relfected across 
% the antenna 
%o 
Et=Ei+Er; 
magt=abs(Et); phast=phase(Et).*(180/pi); 
if loop == 1, 
magihold=magi; magrhold=magr; magthold=magf; 
phashold=phas; phasrhold=phasr; phasthold=phast; 
end 
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% ** END OF MAGNITUDE AND PHASE 





% eer Tre nS rrrrrrr rrr rrrrrrr ort SPITE TEST LET E LCL EL EL LL kL ES hh ke Rieke ieieieieReieieleael 


2G. COMPLEX INTEGRATION OF mag(Et) 


ke 


% de ee HK HT IE KKK TIT IKK KIKI TIKI KIRK IK KKK KIKI KEIR EEE KEKE EKER EN 
1°] 


% 

for m=1:incr, 
for n=1:incr, 
value(n)=Et(m,n)*(d/iner)*(d/incr); 
end 3 

temp1(m)=sum(value); 

end 

field=sum(temp1); 

field=abs(field); 

Total(loop)=20*log10(field); 

end 

% 

Loss(freq)=Total(1)-Total(2); 

end 

%o 

% 

% ** END OF COMPLEX INTEGRATION 


% *™* PLOTTING 


% errr rrr rrr rrrrrrrrrrrrrrrerrre rr eT EPPS LIST TESS Te SL LPL ek eel 


kk 


% errr rrr rrrrrrrrrrerrrrrrerrrerr rE PTS PTE PEEP Sieh tet 2k eee ee aiid 


er=num2str(erx); 
erf=nume2str(ertfx); 
tand=numé2str(tandx); 
tandf=num2str(tandfx); 
h=numestr(h); 
l=num2str(l); 
di=num2str(d1); 
d2=num2str(d2); 

Yo 
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% If range of frequencies were specified plot Loss vs. frequency 
%o 
if count > 1, 
freq=linspace(fstart/1e9, (fstop+1e9)/1e9,count); 
plot(freq,Loss,'*’,freq,Loss) 
xlabel(‘freq (GHz)’) | 
ylabel('Loss (dB)’) 
title(['er= '‘,er,', erf= ‘erf,', tand= ‘tand,', tandf= ',tandf.... 
he Uh,', dt=d2= ',d1,°, l= ‘,I]) 
if pflag == (‘y’), print; end 
end 
% 
% If one frequency was specified mesh the magnitude and phase 
% of DIRECT,REFLECTED and COMBINED rays. 
if count == 1, 
fig=1; 
fs=num2str(f); 
figure(fig); 
meshce(magihold); 
title([‘Magnitude of DIRECT rays wall - f = ‘,fs]); 
xlabel(‘Antenna (incr)'); ylabel(‘Antenna (incr)’); 
zlabel('abs(Ei)’); 
text(.5,1,['er=,er,', erf='erf,', tand='tand,', tandf=',tandf.... 
' he'h,', dt=d2="',d1,', [=',!],'units','normalized’,‘horiz’, center’); 
%e print -depsc plot1_25 
if pflag == ('y'), print; end 
%e 
if ptest ~= 'n’, 
fig=tig+1; 
figure(fig); 
meshc(phashold); 
title((‘Phase of DIRECT rays - f = '‘,fs}); 
xlabel(‘Antenna (incr)'); ylabel(‘Antenna (incr)'); zlabel(‘phase(Ei)’); 


text(.5,1,['er=',er,', erf="erf,’, tand=',tand,’, tandf=',tandf.... 
‘ h='h,', d1=d2=",d1,', l='!],‘units','normalized’,‘horiz','center'’); 
% print -depsc plot2_25 
if pflag == (‘y’), print; end 
end 
% 
fig=fig+1; 
figure(fig); 
meshc(magrhold); 
title([‘Magnitude of REFLECTED rays - f = ',fs]); 
xlabel(‘Antenna (incr)'); ylabel(‘Antenna (incr)'); zlabel(‘abs(Er)’); 
text(.5,1,['er=‘er,', erf=',erf,', tand=',tand,', tandf=',tandf,... 
' h='h,', d1=d2=',d1,', I=',I],'units','‘normalized',‘horiz’,'center'’); 
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% print -depsc plot3_25 
if pflag == (‘y’), print; end 
% 
if potest ~= 'n’, 
fig=fig+1; 
figure(fig); 
meshc(phasrhoid); 
title(('Phase of REFLECTED rays - f = ',fs}); 
xlabel(‘Antenna (incr)'); ylabel(‘Antenna (incr)'); zlabel(‘phase(Er)’); 
text(.5,1,[‘er=",er,', erf=erf,’, tand=',tand,’, tandf=',tandf.... 
' h='h,', d1=d2=',d1,', l='\I],'units','normalized',‘horiz','center’); 
% print -depsc plot4_25 
if pflag == (‘y'), print; end 
end 
% 
if ptest ~= 'n’, 
fig=fig+1; 
figure(fig); 
meshc(phasthold); 
title((‘Phase of COMBINED rays - f = ',fs}); 
xlabel(‘Antenna (incr)'); ylabel(‘Antenna (incr)'); zlabel(‘phase(Et)’); 
text(.5,1,['er=",er,', erf=',erf,’, tand=',tand,’, tandf=',tandf.... 
' h='h,', dt=d2=',d1,', l='\I],'units', ‘normalized’, ‘horiz’, center’); 
% print -depsc plot5_25 
if pflag == (‘y'), print; end 
end 
% 
fig=fig+1; 
figure(fig); 
meshc(magthold); 
Loss=num2str(Loss(freq)); 
title([/Magnitude of COMBINED rays - f = ',fs}) 
text(.5,.9,['Loss w/r free space= ',Loss,'dB'],'units',‘normalized.... 
‘horiz’,'center’); 
text(.5,1,['er=',er,', erf=',erf,', tand=',tand,’, tandf=',tandf.... 
' h='h,', d1=d2=',d1,', l=',I],'units', ‘normalized’, ‘horiz','‘center’); 
xlabel(‘Antenna (incr)'); ylabel(‘Antenna (incr)'); zlabel('abs(Et)’); 
% print -depsc plot6_25 
if pflag == (‘y’), print; end 
end 
% 
% 
% ** END PLOTTING oo 
Of, -anacwvasnanananenenennnoweuenecnaceesnesesnnnnnnnacccnnacacccewecenensannanssenqenseanaenceneeensecese=s 
%o 
% KKKKKKEAKKEKKEKKCREKEKKKKKEEE END PROGRAM KKKKKKKKKKKKKKRKKEKKKKKK KKK KKK 


% were rrerrC CCC CCC SPTCCCLCOOC SSCS SSS SLES LESSEE ET SPSS ETT ET CLC PCT FC CLERC PLPC PARP Ee el 
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APPENDIX C 
SINGLE FREQUENCY PLOTS 
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Figure 51. Magnitude of direct rays at 2.5 GHz 
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Figure 52. Phase of direct rays at 2.5 GHz 
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Figure 53. Magnitude of reflected rays at 2.5 GHz 
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Figure 54. Phase of reflected rays at 2.5 GHz 
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Figure 55. Magnitude of combined rays at 2.5 GHz 
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Figure 56. Phase of combined rays at 2.5 GHz 
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APPENDIX D 
ANTENNA PARAMETER PLOTS 
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Figure 57. Slot antenna impedance measurement - 8 to 12 GHz 
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